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Abstract 

We propose a scalable method on the basis of nth-order coupling operators to construct f- 
dependent phase transformations in the n-qubit modified Deutsch-Jozsa (D-J) quantum algorithm. 
The novel n-qubit entangling transformations are easily implemented via J-couplings between 
neighboring spins. The seven-qubit modified D-J quantum algorithm and seventh-order coupling 
transformations are then experimentally demonstrated with liquid state nuclear magnetic reso- 
nance (NMR) techniques. The method may offer the possibility of creating generally entangled 
states of n qubits and simulating n-body interactions on n-qubit NMR quantum computers. 

PACS numbers: 03.67.Lx, 76.60.-k, 75.10.Jm 
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By using the characteristics of quantum mechanics, quantum computers (QC) are faster 
than their classical counterparts when performing certain computational tasks such as fac- 
torizing a large number |j| , searching unsorted database [2] , and especially when simulating 
quantum systems themselves [3]. Among all quantum algorithms, Deutsch's algorithm QQ 
was not only the first to be proposed, but also the first to be demonstrated p using liquid 
state nuclear magnetic resonance (NMR) techniques 0,0]. An experimental realization of 
Deutsch-Jozsa (D-J) algorithm has not only demonstrated that QCs have an information- 
processing capability much greater than that of their classical counterparts, but also provided 
crucial insight into more impressive quantum algorithms |9j. Based on the Cleve version £|, 
two-, three- or five-qubit D-J algorithm has been implemented with NMR spectroscopy 6|. 
All of the previous D-J algorithms required n+1 qubit quantum systems to realize the n-bit 
D-J algorithm. However, Collins et al. proposed a modified version of the D-J algo- 
rithm, which allowed an implementation of the n-qubit D-J algorithm by using n qubits. Up 
to now, the modified D-J algorithm has been implemented only on three-qubit spin systems 

0. 

Extending qubit number in quantum computation has been a subject of much recent in- 
terest Currently only NMR technique is used to experimentally demonstrate quantum 
information processing with more than four qubits. Such progresses Il2j include seven-qubit 
"cat" -state 131 and Shor's algorithm Il4|, five-qubit D-J algorithm lol, five-qubit order find- 

~i n 

ing algorithm and five-qubit error correction benchmark However, experimentally 
it is a technical challenge to extend to more qubits because of the low signal-to-noise ratio 
(SNR) in NMR experiments and furthermore, the SNR's exponential decrease with the in- 
crease of qubit number In the way to get more qubits, there are two obstacles. First, 
although efficient methods have been developed to refocus the interaction of single-spin and 
two spins 3], it is not easy to control the coherent evolution between two specific spins. 
This is because the coupling network becomes more complex when more qubits are involved. 
Second, it is difficult to find any molecule that has a ppr eciable J-coupling configuration be- 
tween any pair of spins. Fortunately, Madi et al. have used quantum swap gates to 
surmount the second difficulty 

In this letter, we propose a simple method to construct n-qubit entangling transforma- 
tions on the basis of nth-order coupling operators. The method we present is general and 
scalable. The key to the method is only to use the J-couplings between neighbouring spins to 
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implement nth-order coupling transformations. With this method, we can construct n-qubit 
entangling transformations for balanced functions in the modified D-J algorithm. And then 
we experimentally demonstrate the seven-qubit modified D-J algorithm and seventh-order 
coupling transformations using NMR. 

The modified D-J algorithm. — A D-J problem 0, is: given Boolean function f(x): 
{0, l} n — > {0, 1}, how many evaluations are necessary to determine f(x) being constant or 
balanced — constant functions f(x) always give the same output (all or all 1) for all input 
values, but balanced functions f(x) have an equal number of outputs as Is. In a modified 
version of the D-J algorithm [10], the key to the scheme is to perform an f-controlled gate 



whose transformation on the basis elements \x >= |x„_i • ■ ■ x > [xj G (0, 1)] is defined 



id 

as 




1( 



U f \x>:={~l) f{x) \x> . (1) 

The unitary transformation Uf is an f-dependent phase operation on 2 n eigenstates of an 
n-qubit quantum computer. For example, the number of functions for the seven-bit D-J 
problem is (g 2 4 8 )-|-2, which is a large number. In general, it is difficult to express the unitary 
transformation Uf as the product of operators allowed by the Ising-type Hamiltonian of 
spin-1/2 systems. Fortunately, we can classify the Uf s as seven representative types of 
unitary transformations namely, non-entangling, and from two-qubit entangling to seven- 
qubit entangling. 

Constructing n-qubit f-dependent phase transformations via nth-order coupling 
operators. — We consider 9 explicit transformations, as shown in table 1. The transforma- 
tions (Uf4 — Uf 9 ) in table 1 correspond to the two-qubit to seven-qubit entangling operators, 
respectively. From the transformations (UfA — Ufg), we can easily extend the entangling 
transformations to an n-qubit case. For example, an f-dependent phase transformation for 
any balanced function in the n-bit D-J algorithm can read: 

= \ (E {1) <8> • • • <g> E {n ~ 2) <g> ® E(n) + E(l) ® • • • ® £ ( ™~ 2) ® ® ^ n) ) 

A ® . . . <g, a f- l) <g> £7< n > - (jW ® • • • ® a^) , (2) 
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where o z is the Pauli matrix, E is a 2x2 unit matrix, and superscripts label the qubits in- 
volved. The transformation U n -ut-f is diagonal in the 2 n eigenbasis vectors. It is interesting 
to note that the transformation U n _bit-f cannot be decomposed as direct products of single- 
spin operators and is hence an entangling operation. The n-qubit entangling transformation 
U n -bit-f in the modified D-J algorithm can be constructed on the basis of second, nth, and 
(n-l)th-order coupling operators. In order to realize the n-qubit entangling transformation 
U n -bit-f experimentally, U n -bu~f can be decomposed with the product operator basis set: 
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U n - b it-f = exp[in\ |exp 

E {1) ® • ■ • ® E {n ~ 2) <g> exp 
E {1) ® • ■ • ® E [n - 2) ® exp 



(1) . . . 



E {n) 
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Following a procedure similar to that given in Ref. 



we perform the nth-order cou- 



pling operation exp 



l 2 a 



(1) 



(n) 
0"z 



in terms of the operators allowed by the Hamiltonian 
of spin- 1/2 systems, i.e., with radio frequency pulses and scalar J-couplings between two 
neighboring spins. For example, a quantum circuit for performing a seven-qubit entangling 
transformation (£7/9 ) is demonstrated in Fig. 1(a). 

Experimental realization. — A quantum circuit shown in Fig. 1(b) can be used to im- 
plement the seven-qubit modified D-J algorithm. This circuit begins with an initial state 



(|0000000>). In fact, however, the statistical mixture of pure states can be used |6l llOL |21| 



as input states, since the D-J algorithm can accept the thermal equilibrium state as an input 



2l|. There are two important operations in Fig. 1(b), namely pseudo-Hadamard and 



f-dependant phase gates, h and h~ x are pseudo-Hadamard gates, which can be implemented 
by 90 degree pulses along ±y axes, respectively. The key to this approach is to realize 
seven-qubit f-dependant phase transformations shown in table 1 and Fig. 1(a). 

In order to read out results, it is necessary to apply a selective (vr/2)^ pulse to the 
corresponding spin. According to phases (absorption or emission) of the 13 C and 1 H NMR 
signals, NMR spectra clearly indicate where the system is. Therefore we can determine the 
constant or balanced nature of the function (Fig. 1(b)). Fortunately, we can cancel the final 
pseudo-Hadmard transformations by the (vr/2) J/ read-out pulse used in NMR experiments. In 
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order to improve experimental results, we have used the first-order phase correction methods 
2^|, and acquired NMR spectra with 96 scans. 

We have selected a seven-qubit spin system [U- 13 C 4 -labeled crotonic acid ( 13 C 1 H| 
i3 C 2 H i = i3 C 3 H 2 i3 C 4 02 H, Cambridge Isotope Laboratories Inc. Cat. No. CLM-6118), 
where right superscripts label the qubits] |l3[. The sample of 20mg crotonic acid was dis- 
solved in deuterated acetone, degassed and flame-sealed in a standard 5mm NMR test tube, 
the coupling constants of which can be found in Ref. We began the experiment from 

a thermal equilibrium state. We have used four carbons and three protons as the seven 
qubits. We carried out the experiment with a Varian INOVA 600 NMR spectrometer. All 
NMR experiments were conducted at 25 °C. The shape of the soft pulse was Gaussian. 
The J-coupling interactions in NMR spin systems have been used to implement entangling 
transformations required for the seven-qubit modified D-J algorithm. 

Results and discussion. — Typical experimental results are shown in figures 2 and 3. The 
experimental spectra in figure 2 correspond to the constant function (E//i), which is directly 
obtained after applying a 90 degree pulse to the thermal equilibrium state. The spectra 
in Fig. 2 served as reference spectra. The phase of the reference spectra was adjusted so 
that signals from all the seven spins appear in absorption (positive phase). The experimental 
spectra in Fig. 3 correspond to the balanced functions (Z7/2, Uf±, Ufs, and U/g), respectively. 
Compared with the reference spectra in Fig. 2, the spectra in Fig. 3 show that there is at 
least one line with a tt phase difference (in emission). By determining the relative phase of 
the signals from the seven spins, we can determine that the functions corresponding to Fig. 
3 are balanced. It should be noticed that only a single function call is used. 

Fig. 2 shows that the spectra of C 1 contain 128 peaks. This indicates that spin C 1 
interacts with 7 neighbors. This results from the J-couplings of the three other 13 C, the two 
protons (H 1 , H 2 ) and the three methyl protons (H 3 ). It should be noted that the J-coupling 
effect of the methyl protons on the spin C 1 is the same as that of two protons. 

As for the thermal equilibrium state used as an input state in our experiments, the 
investigations 0, Q, 0| have shown that the D-J algorithm can be implemented by starting 
with thermal rather than pure initial states. The reason is that the thermal equilibrium 
state has similar effect as the effective pure state. 

There are imperfect phases in the experimental results. The signals do not exhibit pure 
absorption or pure emission lineshapes. The phase errors mainly come from (a) the im- 
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perfection of the selective pulses, (b) the inaccuracy of the 90 and 180 degree pulses, (c) 
inaccurate refocusing of chemical shifts during the J-coupling delays, and (d) J-coupling 
evolution during long selective pulses. 

It should be noted that n-qubit entangling transformations in the modified D-J algo- 
rithm can be performed with radio frequency pulses and scalar J-couplings between two 
neighboring spins. The scheme may be used to create generally entangled states with n 
qubits on the basis of the interactions between neighboring spins. This is very important 
for demonstrating the role of quantum entanglement in quantum information processing. 

We used seventh-order coupling operators to construct a seven-qubit entangling trans- 
formation. We have successfully realized the seven-qubit entangling transformation in the 
context of the modified D-J algorithm. Therefore we have experimentally performed the 
seventh-order coupling operator. In fact, there is no seventh-order coupling in liquid state 
NMR systems. Our experimental results show that seven-body interactions can be simulated 
with such an experimental method. This will be useful in simulating many-body interactions 
in other quantum systems. 

The NMR experimental realization of multi-qubit quantum algorithms is a difficult propo- 
sition. In a multi-qubit NMR QC, it is necessary to suppress the interactions coming from 
the other qubits when realizing the quantum logic gate between two given spins. Our exper- 
imental results have demonstrated that the control of multi-qubit's evolution is successful. 

Conclusions. — On the basis of constructing f-dependent phase transformations in the 
modified D-J quantum algorithm with nth-order coupling operators, we have experimentally 
tested the modified D-J algorithm for seven qubits and performed seventh-order coupling 
transformations. This may open a way to experimentally simulate n-body interactions and 
to realize n-qubit entangling transformations on NMR quantum computers. 
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Table 1. The transformations U / and their defintions. E^ 3 ' '"'^ and a{ A ' B ' C '" mean 
E( A ) <g> ® ® . . and cr£ A) <g> a[ B) <g> af ] <g> ■ ■ •, respectively. C i (i=l,2,3, and 4) and 
FT (j=l,2, and 3) are carbon and hydrogen nuclei used in experiments. 



U f Defintions 



Constant: 

U E {C\C 2 ,C Z ,C\H\H 2 ,H 3 ) 

Balanced: 

U f2 a f\C*) ^ E{C ^C\H\H\H*) 

jr (C\C\C\C\H\H\H*) 

U f4 \ {E^^) + af 3) ® E^) + ® aP - vf^) ® <rP ® i^'W 3 ) 
| (i^) ® erP ® E^ 3 ) + ® af^ HZ) + 2 ' C1) ® E^ 3 ) - of 2 > C1 ' H3) ) 

I ^(C 1 ^) a f 3 ) + ^.C) a f 3 ^ 4 ) j E {H\H*,HS) 

+ \ (af^ C3) ® E^ - of ® £(^ 2 ," 3 ) 

Uf7 +1 (<rf 1 ' C2 ' C3 ' C4) ® E^ 1 ) - a f 1 ^ 2 ^ 3 ^ 4 ^ 1 )) £(" 2 ," 3 ) 



f7 



/5 



C 2 ) 
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Figure captions 

FIG. 1. (a) A quantum network for realizing the seven-qubit entangling transformation t//g. 
Refocusing and decoupling pulses are not included. Vertical bars denote J-coupling gates, 
which are given by the unitary transformation exp{— ia^a^n/A} on spins k and 1. (b) A 
quantum circuit for implementing the seven-qubit modified D-J algorithm. A function f is 
constant if and only if the output state for every qubit is in all the |0> state, otherwise a 
function is balanced. 

FIG. 2. The experimental spectra corresponding to the transformation (£//i) for the con- 
stant function fl from seven spins. Horizontal axes, relative frequency (ppm); vertical axes, 
intensity (arbitrary units). 

FIG. 3. Same as Fig. 2, but (a), (b), (c), and (d) corresponding to the balanced function 
transformations Uf2, Uf4, C//5, and C//9, respectively. 
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